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ABSTRACT 

An overview of the results of the Near-IR SO galaxy Survey (NIRSOS) is presented. 
NIRSOS is a magnitude (mB^12.5 mag) and inclination (< 65°) limited sample of 
~ 200 nearby galaxies, mainly SOs, but include also Sa and E galaxies. It uses deep 
ii's-band images, typically reaching a surface brightness of 23.5 mag arcsec"^. De- 
tailed visual and photometric classifications were made, for the first time coding also 
the lenses in a systematic manner. The main analysis methods include 2D multi- 
component decomposition approach, and Fourier analysis of the non-axisymmetric 
structures. As a comparison sample, a similar sized spiral galaxy sample with similar 
image quality was used. The main emphasis were to study whether the SOs are for- 
mer spirals in which star formation has been ceased, and also, how robust are bars in 
galaxies. Based on our analysis the Hubble sequence was revisited: following the early 
idea by van den Bergh we suggested that the SOs are spread throughout the Hubble 
sequence in parallel tuning forks as spirals (SOa, SOb, SOc etc.). This is evidenced by 
our improved bulge-to-total (B/T) flux ratios in the SOs, covering the full B/T range, 
reaching small values typical to late-type spirals. The properties of bulges and disks 
in SOs were found to be similar to those in spirals and also, the masses and scale 
parameters of the bulges and disks to be coupled. It was estimated that the spiral 
bulges brighter than -20 mag in K-hand are massive enough to be converted into the 
bulges of SOs merely by star formation. Bars were found to be fairly robust both in 
SOs and spirals, but inspite of that bars might evolve significantly within the Hubble 
sequence. 
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1 INTRODUCTION 

We review the main results of the Near-IR SO galaxy Survey (NIRSOS) obtained so far0. NIRSOS is a magnitude (m_B^12.5 
mag) and inclination (less than 65°) limited sample of ^ 200 nearby galaxies, mainly SOs, but include also Sa spirals and 
25 late-type ellipticals. Late-type ellipticals were included for not to miss any potentially misclassified SOs. The observations 
were done in the A"s-band, carried out using 3-4 meter sized ground-based telescopes with sub-arcsecond pixel resolution. 
The images are deep, typically reaching a surface brightnesses of 23.5 mag arcsec"'^ in azimuthally averaged profiles (~ 2 
mag deeper than the 2MASS images), thus allowing the detection of the faint outer disks in SOs. Our main emphasis was to 
address possible secular evolutionary processes in galaxies by comparing the photometric properties of SOs and spirals, based 
on similarly selected samples, with similar image quality. 
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In the early classification by Hubble [37], the SOs were an enigmatic group of galax;ies between the ellipticals and early- 
type spirals, and since then they have appeared to be important in any galaxy evolutionary model. SOs are suggested to 
evolve from spirals, cither by internal secular processes in isolation, or by gas stripping mechanisms, followed by quiescent 
star formation ([46]): in cluster environment by ram pressure stripping ([36]) or by galaxy harassment ([57]), and in galaxy 
groups by tidal encounters (e.g. [38]). Alternatively, SOs have been suggested to form by galaxy mergers in a similar manner 
as elliptical galaxies ([68], [6]), or have accreted most of their mass in minor mergers ([3], [41]). 

Originally, the stripping scenario was adopted because SOs were found to appear mainly in galaxy clusters ([26]), and 
also because they were kinematically more related to spirals than to ellipticals ([27]). However, the large bulge fractions found 
by Simien & de Vaucoulcurs ([67]) in SOs has worked as a counterargument to this scenario, because spiral galaxies do not 
have sufficient amount of interstellar gas to build such large bulges ([45]). Some SOs are also found to have chemically and 
kinematically decoupled cores ([1]), which are generally interpreted in terms of sinking of small gas poor satellites into the 
main galaxies. Likewise, the merger approach has problems, at least if interpreted as the dominant mechanism for producing 
the SOs. For example, modern semi-analytical merger models (e.g. [41]) successfully produce the large bulge masses by Simien 
& de Vaucouleurs [67], but fail to recover the observed trend in a/Fe elements with stellar velocity dispersion ([69], [40]). 
They have also difficulties to explain the observed tight luminosity-size relation for the SOs ([60]). 

The morphological division between E and SO galaxies was obscured by the detection of boxy and disky ellipticals by 
Bender ([12]; see also [44]). It was also found that both the galaxy luminosity and the degree of rotational support [a/Vmax) 
correlate with the isophotal shape of the elliptical galaxies ([13]). Among other things, these findings led King ([42]) and 
Djorgovski ([25]) to announce that the Hubble sequence is breaking down and should be replaced by a more physical approach, 
based on measured physical parameters, related to galaxy kinematics. Recently it was suggested by Cappellari et al. [22] that 
such a parameter could be Xrc, a proxy for the specific angular momentum in galaxies. Based on Are they showed that 
most SOs, and also even 2/3 of the elliptical galaxies in the nearby Universe are rotationally supported. It was also shown 
by Emsellem et al. [29] that this kinematic parameter has no clear correlation with the boxy /disky isophotal shapes. These 
findings arc suggested to put into a now light many of the previous results, including the morphological classification of the 
early- type galaxies ([22]). NIRSOS has a large overlap with the sample by Cappellari et al. [22] (used in [29]), thus allowing 
to compare the kinematic classification with the new morphological classification of SOs made by Laurikainen et al. ([54]). 

In the following an overview of the main results of NIRSOS is given, discussed in conjunction of some recent studies of SOs 
(more complete references can be found in the original papers). As a comparison sample, the Ohio State University Bright 
Spiral Galaxy Survey (OSUBSGS; [31]) has been used, for which galaxies we have carried out similar analysis as for NIRSOS. 
Our main analysis methods include detailed visual and photometric classification of galaxies, 2D multi-component structural 
decompositions, Fourier analysis for calculating the properties of bars, and isophotal analysis. In this paper we also analyze 
a sub-sample of NIRSOS, which have kinematic information given by Emsellem et al. ([29]). A picture is outlined supporting 
the view that SOs in general are evolved from spiral progenitors. They are suggested to be spread throughout the Hubble 
sequence forming SOa, SOb, SOc types, in a similar manner as spirals. Our scenario, originally outlined in Laurikainen et al. 
[53], is similar to that suggested recently by Kormendy & Bender [88]. 



2 NIRSOS ATLAS 

NIRSOS Atlas ([54]) presents the flux calibrated images, shown in a logarithmic scale. Also shown are the dimensions of the 
identified structures, as well as the radial profiles of the isophotal analysis results (ellipticity, position angle, and the parameter 
b4, which measures the isophote's deviations from perfect ellipticity). An example of the layout of the atlas images is shown 
in Figure 1. 

2.1 Morphological classification and how it compares with kinematic clcissification 

The classification is based on de Vaucouleurs' revised Hubble-Sandage system ([73], see also [19], [21]), which includes the 

stage (S0~, S0°, SO"*", Sa), family (SA, SAB, SB), variety (r, rs, s), outer-ring or pseudo-ring (R, R'), possible spindle shape, 
and the presence of peculiarity. The ring classification follows Buta & Crocker ([15]) and Buta ([16]). What is new in our 
classification is that lenses (nuclear, inner, outer) are, for the first time, systematically coded for a significant sample of SOs. 
A new lens type, 'barlens', is also introduced: it is a lens-like structure embedded in a bar, with a fairly sharply declining 
outer surface brightness profile - this component has been often erroneously mixed with the bulge. Barlenses are shorter than 
the main bars, and are typically identified inside strong bars. Due to the sub-arcsecond pixel resolution it was also possible 
to classify the central structures like nuclear bars, rings and lenses in a systematic manner. 

Besides the above discussed visual classification, we also made a photometric classification where faint structures are 
identified even if they were not directly visible in the images. Faint structures, overshadowed by bulges, were detected after 
subtracting the bulge model in the decomposition, and/or by making unsharp masking. Lenses were detected by inspecting 
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the surface brightness profiles, in which they generally appear as nearly exponential subsections, however with more shallow 
profiles than the main disk component. Photometric classification turned out to be critical for separating SOs from the elliptical 
galaxies at the E/SO interface. 

NIRSOS has 66 galaxies in common with the SDAtlas, which is a volume-limited sample of 260 early-type galaxies, with 
the kinematic Ajje parameter available ([29]). Ajje is a proxy for the specific angular momentum in galaxies, and can be used 
as a measure for the rate of the kinematic support, thus dividing the galaxies to fast and slow rotators. It is measured within 
one effective radius of galaxy brightness thus covering well the bulge region. Of the 66 galaxies in common, 9 were ellipticals, 
and 56 were SO-SO/a galaxies in the Third Reference Catalog of Bright Galaxies ([73], RC3). However, in our classification 
only 2 out of the 9 galaxies were truly ellipticals. Kinematically, one of the two ellipticals is a slow rotator, whereas the other 
is a fast rotator showing an inner disk. The complete NIRSOS has 26 galaxies classified as ellipticals in RC3, of which only 3 
remained ellipticals in our photometric classification. This comparison shows that generally the Hubble stage, based on the 
morphological and kinematic classification is consistent with each other (notice that NIRSOS does not include highly inclined 
galaxies). However, based on Xrb not all galaxies with prominent outer disks are fast rotators. Slow rotators include NGC 
4552, NGC 5846, NGC 5631 and NGC 6703, in which the SO nature is confirmed by the detection of a lens (NGC 5846 even 
has multiple lenses). Slow rotators comprise also NGC 4406 and NGC 4472, which have exponential outer profiles, used as 
evidence of a disk in our photometric classification. In conclusion, this means that Xrs alone cannot be used to distinguish 
the Hubble stage among the SOs and the elliptical galaxies. 

3 PROPERTIES OF BARS 

3.1 Morphology of bars in SOs 

Bars are the main characteristic common in SOs and spirals. They can be regular elongated structures (NGC 4608), have 

condensations at the two ends of the bar (ansae; e.g. NGC 2983), or show X-shaped structure in the inner part of the bar 
(e.g. IC 5240) (Fig. 2). Besides ansae, NGC 2983 and NGC 4608 have also barlenses. Both structures are illustrated in Figure 
2 by subtracting the underlying bar+bulge model from the image of NGC 2983: the residual image shows a nearly spherical 
barlens, and two blops (ansae) that form part of the bar. 

It was shown by Laurikainen et al. ([51]) that ansae in bars are characteristic to SOs, appearing less frequently in spiral 
galaxies (40% and 12%, respectively). This was shown simultaneously also by Martinez- Valpuesta, Knapen & Buta ([55]), who 
additionally suggested that strong bars (B) have ansae more frequently than weak bars (AB). Both studies used a sub-sample 
of NIRSOS, and in [55] also the de Vaucouleurs' atlas by Buta et al. ([19]). By repeating the statistics for the complete NIRSOS 
sample in this study, using the classifications of the NIRSOS Atlas, we find ansae in 30% of the SOs, which is still clearly higher 
than in spirals. However, now the fraction of ansae is the same for strong (B) and weak (AB) bars. In Laurikainen et al. [52] 
the large fraction of ansae bars among the SOs was discussed as evidence of bar evolution in the Hubble sequence. 

In NIRSOS Atlas X-shaped bars appear in 9 galaxies, which is 8% of all barred galaxies ([54]). This is important, because 
all the galaxies in NIRSOS have inclinations less than 65°, although the X-shaped bars are generally discussed in the context 
of nearly edge-on galaxies, suggested to be formed by vertical thickening due to buckling effects ([23], [79]). It was shown by 
Athanassoula & Beaton [7] that the vertical thickening might still be recognized at inclinations of 77.5°, whereas our study 
shows that the X-shaped bar morphology can be visible even in almost face-on galaxies. A nice example of such a bar appears 
in IC 5240 (Fig. 2), which has an inclination of 49" (taken from [54]). Evidently, bar buckling alone, leading to a vertical 
thickening of a bar, is not sufficient to explain such bar morphologies. 

Multiple bars appear in ~ 20% of the SOs [52], which is a similar number as found previously also for spirals [75, 76]. It 
is worth noticing that SO/a galaxies have a larger number of multiple bars [52]. A prototypical case where two bars appear 
nearly perpendicular to each other is shown in Figure. 3, in which both the nuclear and the main bar are surrounded by 
lenses. Some galaxies, like NGC 2681, have even three bars. 

3.2 Fourier properties of bars: SOs vs. spirals 

The properties of bars in SOs (using a sub-sample of NIRSOS) and spirals (using OSUBSGS) has been studied by Laurikainen 
et al. [51] using a Fourier method. For the complete NIRSOS they will be presented in a forthcoming paper. Tangential forces 
induced by the non-axisymmetric structures, mainly bars, were measured using a bar torque approach ([24]), applying the 
polar method by Salo et al. ([62], [64]; see also [47]). In this method the surface brightness at each radius is decomposed into 
azimuthal Fourier components. Using the even Fourier components (typically up to m=10), the gravitational potential in the 
equatorial plane of the galaxy is obtained by FFT over azimuth, combined with direct summation over radial and vertical 
coordinates. An important advantage of the polar method is the suppression of spurious force maxima which may arise in the 
direct 3D Cartesian FFT integration in the noisy outer disks. The calculated Fourier amplitude profiles also provide another 
useful tool for characterizing bars (see below). 



4 E-mail:eija. laurikainen@oulu.fi 



As a measure of bar strength, we use Qg = max{FT/ < Fr >), the maximum of tangential force amplitude relative 
to the mean axisymmetric radial force, evaluated at the region of the bar. The mass-to-light ratio {M/L) is assumed to be 
constant, and the vortical scale height of the disk (and bar) is estimated from the exponential scale length, using a Hubble typo 
dependent mean ratio. Also, the different 3D density distribution of the bulge is corrected, based on bulge models obtained 
from decompositions. The effect of including dark halo force field was also investigated, but its influence on Qg appeared to 
be insignificant at the bar region ([63], [17]). These calculations also give a proxy for the bar length, rQg, which is the radius 
where the max;imum tangential force Qg occurs. Bar lengths were estimated also from the phases of the A2 Fourier amplitude 
by assuming that it is maintained nearly constant in the bar region (a correlation between this bar length and rQg was 
shown by [48]). The maximum of m=2 Fourier amplitude, A2, was used as an estimate of the relative brightness of the bar. 
These properties are shown in Figure 5 in Laurikainen et al. [51] as a function of Hubble: bars grow in length and in relative 
brightness (^12) towards the early- type galaxies, but for Qg the trend is opposite. Notice that although the bar ellipticity 
(shown in the same figure) correlates with Qg ([91]), it has no systematic correlation with the Hubble type. In Laurikainen 
et el. [48 ] the tendency of weakening bar strengths {Qg) towards the early- type galaxies was explained by a dilution effect 
due to the more massive bulges in the early-type galaxies (e.g. the average Qg parameter may decrease even if the average A2 
amplitude increases, since the bulge contribution to radial force becomes more important toward earlier types). There exist a 
correlation also between Qg and A2, but for the above reason the correlations are different for the early and late- type galaxies 
(see Fig. 8 in [48]). The obtained tendency for bar lengths was originally shown by Elmegreen & Elmegreen [28]. 

For a sub-sample of 26 barred galaxies in NIRSOS, the radial A2 profiles were fitted by single (SG) and double (DG) 
Gaussian functions by Buta et al. [18]. It appeared that 65% of the bars in SO-SO/a galax;ies have single Gaussian profiles, 
whereas 35% are best fitted by two Gaussian functions. Typical examples of such profiles are shown in Figure 4. The galaxies 
with DG bars typically have also significant higher Fourier modes in the bar region (m=4, 6, 8), in addition to m=2 ([51]). 
It was discussed by Buta et al. that the DG-profiles are similar to those predicted by the simulation models (e.g., [5]) in 
which the bar transfers a large amount of angular momentum to the halo. An attempt to associate the DG-profiles to specific 
morphological structures was made by Laurikainen et al. [51] who suggested that the fat or double-peaked Gaussian amplitude 
profiles are due to two bar components, a long and narrow bar, and a shorter component in the inner parts of the bar (or 
an inner oval). DG bars were found to be more prominent, not only in terms of Qg, but also in A2 and bar length [51]. In 
Laurikainen et al. [54] those inner bar components were associated mainly with barlcnsos (though some of them can be ovals), 
which are found to appear in 30% of barred SO-SO/a galaxies in NIRSOS. A good example is NGC 4314 (Fig. 1), in which 
the barlens is the fat elongated structure inside the bar. Erwin et al. [30] have discussed two SOs, in which a superposition of 
a classical and a pseudo-bulge was suggested. These galaxies are NGC 2787 and NGC 3945, which form part of NIRSOS. In 
both galaxies the component interpreted as a pseudo-bulge by Erwin et al., is called as a fat inner bar component in [51], and 
more recently defined as a barlens by us [54]. 

3.3 Do SOs have the bar strengths expected for systems not accreting any gas? 

The above question was recently made by Buta et al. [86] with the main emphasis to test the hypothesis by Bournaud & 
Combes [14], in which multiple bar episodes are expected in the Hubble time. In this scenario bars form and evolve in galaxies 
when they have gas, and the evolution stops when the gas in used in star formation. These stars are then transferred into 
the bulge, for example by bars or spiral arms in the central regions of the galaxies. When the central mass concentration 
formed by star formation becomes very high, the bar will be destroyed. Therefore, if bar strength varies over time the relative 
frequency of galaxies in each Qg bin tells us the relative amount of time a galaxy spends in a certain bar state (strong, weak, 
non-barred). For spirals this was first tested by Block et al. [80]. They suggested that galaxies might have doubled their mass 
in 10^° years (see Fig. 5, left panel), evidenced by the extended tail towards strong bars, and the lack of weak bars, which 
features are predicted in the strong gas accretion models by Bournaud & Combes [14]. 

This test was later repeated by Buta, Laurikainen and Salo [87] for the same galaxy sample, but using the refined bar 
torque method described in the previous section (we also discussed why the obtained Qg distribution was different from that 
by Block et al. [80]). In Buta et al. [17] the bar and spiral fluxes were additionally separated firom each other. Although the 
correction affected Qg in a few individual cases having very strong spiral arms, it barely affected the Qg distribution. The 
refined Qg distribution ([87], [17]) has a large number of weak bars lacking from that obtained by Block et al. Likewise, it 
has a slightly smaller number of very strong bars. In fact, the obtained Qg distribution (see Fig. 8a in [17]) largely resembles 
the non-accretion model by Bournaud & Combes shown in Figure 5 (left panel) . thus supporting the view that bars in spirals 
are fairly robust. In Buta et al. [86] the Qg distribution for NIRSOS was calculated. Most importantly, a clear difference was 
found between SOs and early-type spirals (see Fig. 5, right panel). This was suggested to support the view according to which 
SOs have not accreted gas for a long time, evidenced by the lack of the extended tail, and the existence of a large number of 
weak bars. As discussed above (see Section 3.2) the smaller number of strong bars among the SOs can be due to a dilution 
effect caused by the more massive bulges and thicker disks in SOs. However, it was also discussed by Buta et al. [86] that this 
cannot produce all of the difference in Qg between the SOs and early-type spirals: although spirals have a larger number of 
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strong bars than SOs, still there are no galaxies having Qg > 0.5. The conclusion in this study was that if SOs are stripped 
spirals, the weaker bars in SOs could indicate that bar evolution continues to proceed even after gas depletion in galaxies. 



4 LENSES 

Lenses appear as flat disk components with rather sharp outer edges ([43]). However, not all lenses are directly visible in 
the images. In the NIRSOS Atlas ([54]) lenses were generally detected as exponential subsections in the surface brightness 
profiles. NGC 524 (Fig. 6) shows all the main lens types, nuclear (nl), inner (1), and outer lens (L). When the outer lens is 
very prominent compared to the underlying disk, as in NGC 1533, it manifests as a broad bump in the surface brightness 
profile, in this case having also some characteristics of a ring (RL) . 

For clarity the examples shown of the different lens types are for non-barred galaxies. However, lenses appear both in 
barred (61%) and non-barred (38%) SOs, based on the classification in NIRSOS Atlas [54]. In Laurikainen et al. [52] even a 
larger fraction of lenses was found, but it was based on a sub-sampic of NIRSOS. In barred galaxies nuclear (nl) and inner 
(1) lenses typically end up to the radius of the nuclear and the main bar, respectively, relating them to resonances of the 
rotating bar (see NGC 1543 in Fig. 3). However, not all lenses are related to resonances. For example, series of lenses in some 
non-barred SOs appear, like in NGC 1411, which cannot be immediately understood in the framework of the resonance theory. 
Lenses can also be relics of significant star formation in the spiral arms as suggested in the NIRSOS Atlas ([54]). Originally the 
idea is from Bosma [8], who also confined the lens formation to the epoch of galaxy formation. Lenses are also suggested to 
form by disk instability, in a similar manner as bars [77] . It was further suggested by Kormendy [43] , that bars may gradually 
dissolve into lenses. In fact there are many results in NIRSOS which are consistent with this scenario: (a) lenses in barred SOs 
often end at the bar radius, (b) SOs were found to have a smaller bar fraction, and a larger fraction of lenses than spirals 
([52]). Also, (c) dissolution of bars would explain the large number of lenses in non-barred SOs in a natural manner. Using 
the ellipticity of a bar, a smaller bar fraction in SOs, compared to that in spirals, was found also by Aguerri, Mendez-Abrey 
& Corsini [4]. Most probably lenses have multiple origins, and in order to better understand their nature detailed analysis of 
their dimensions and physical properties needs to be performed. A forthcoming NIRSOS paper will focus on that. 

It is worth noticing that multiple lenses appear even in 25% of the SOs in the NIRSOS Atlas, including barred and non- 
barred galaxies [54], which needs to be understood in the formation and evolution of galaxies. For example, if a large fraction 
of the mass in the SOs was accreted by minor mergers, it needs to be understood how the multiple lenses can survive through 
such processes. 



5 2D MULTI-COMPONENT DECOMPOSITIONS FOR NIRSOS 

5.1 The multi-component approach 

A 2D multi-component code, BDBAR (written by Salo, and described in [48], [49]), was used for decomposing the light 
distributions of the i^'s-band images into bulges, disks, bars, ovals and lenses. This multi-component approach turned out to 
be important, not only for barred galaxies, but also for galaxies with prominent lenses. Using artificial images this was tested 
by Laurikainen et al. [50]. Figure 7a shows the surface brightness profile of a synthetic image with a bulge and a disk, with 
random noise added, whereas 7b and 7c show the same image after adding a small bar on top of that. Making a simple bulge- 
disk decomposition for the barred synthetic image (in the middle) overestimates B/T (=0.36), due to erroneous assignment 
of the bar flux to the bulge, whereas the bulge-disk-bar decomposition (right) recovers the correct B/T value {B/T=0.27). 
The residual image also shows a bar in the simple bulge-disk model, but not in the bulge-disk-bar decomposition. 

A test for the observed NIRSOS images was made by Laurikainen, Salo & Buta [49], collected to Table 1, where ID (using 
azimuthally averaged profiles) and 2D decompositions are also compared. It appears that simple bulge-disk decompositions 
give a similar mean i3/T- ratio, independent of whether ID or 2D fitting is used, whereas the three-component approach gives 
significantly lower B/T. The value < B/Tk >=0.55 in the bulge-disk decomposition is very similar to < B/Tb >=0.57 
as obtained by Simien & de Vaucouleurs [67]. We also estimated that the different wavelengths used does not cause this 
difference. Adding even more components (=" final" model in the Table), like nuclear bars, further lowers the B/T, but the 
change is not as dramatic as between the 2 and 3 component models. The Sersic index is also smaller in the bulge-disk-bar 
decompositions, but even in the simplest models the mean value is not as large as 4, as often produced by merger simulations. 

5.2 Low B/T values obtained both for SOs and spirals 

The 2D multi-component decompositions for the complete NIRSOS sample are discussed by Laurikainen et al. [53], where 
they are also compared with a similar sized sample of OSUBSGS spirals, using the same decomposition approach (see [48], 
[51]). Internal dust correction was applied to all galaxies, in a similar manner as in Graham & Worley [35]. For spirals these 
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Table 1. Mean _B/T-ratio and Sersic n parameter of the bulge as estimated using the ID and 2D-decompositions, applied for the A's-band 
images. This table was originally published in Laurikainen et al. [49] . 



method < B/T > < rij^.^ig^ > 

2D (final) 0.25 ± 0.03 2.1 ± 0.1 

2D (bulge/disk/bar) 0.30 ± 0.03 2.1 ± 0.1 

2D (bulge/disk) 0.55 ± 0.06 2.6 ± 0.1 

ID 0.48 ± 0.04 2.7 ± 0.4 



results are in good agreement with the bulge-disk decompositions made for non-barred galaxies ([78,84]), and with other multi- 
component decompositions ([81], [82], [83]). It is encouraging that Weinzirl et al. [81], applying bulge-disk-bar decompositions, 
obtained very similar low B/T ratios as obtained previously by Laurikainen et al. ([48], [51]) for the same sample of spirals. 
In both studies M/L ratio was assumed to be radially constant. Graham & Worley [35] also obtained fairly low B/T values, 
but their study uses heterogeneous data from the literature, including simple and multi-component decompositions. [f| 

The dust-corrected B/T histograms for the morphological type bins are shown in Figure 8: though finding an overlap 
in B/T among the Hubble types is not new, this figure clearly demonstrates how B/T in SOs extends to the region of the 
late-type spirals, thus covering the full B/T range to near zero values. To our knowledge this is the first observational evidence 
in support of the scenario suggested by van den Bergh [72], according to which the SOs might form a sequence of SOa, SOb 
and SOc galaxies in a similar manner as spirals. 



5.3 Low B/T in conjunction with semi-analytic models (e.g. Khochfar et al. [41]) 

Formation of early- type galaxies (ETC) has been recently modeled by Khochfar et al. [41] by semi-analytical models (SAM), 
taking into account the new 3DAtlas kinematic observations by Emsellem et al. [29], defining the fast and slow rotators. In 
their models B/T is a critical boundary condition to define the ETCs. The range 0.5 < B/T < 0.9 is used for 80% of the 
fast rotators, and 20% of them have B/T <0.5, for which galaxies also the gas fraction was used as a selection criterion. The 
galaxies with low B/T are similar to late-type spirals, called 'red spirals', except that they have barely no gas, and therefore 
have no spiral arms. In their models the majority of fast rotators have accreted their bulge mass in minor mergers, whereas 
the 'red spirals' were formed by starvation in more dense galaxy environments. 

The above models seem to explain well the distribution of Xrc, a proxy of specific baryonic angular momentum content 
in galaxies. The difference between fast and slow rotators is largely based on the stellar disc fractions in galaxies. This is 
consistent with our morphological analysis of NIRSOS in which, besides SOs, also most of the galaxies classified as late-type 
ellipticals in RC3 show disk structures. However, in many SOs also small central components appear in the regions covered by 
the kinematic observations. This complicates the fast/slow interpretation, because nuclear bars, lenses and rings might have 
different velocity dispersions. The richness of the central structures in SOs might partly explain why B/T does not correlate 
with Afle, which is shown in Figure 9 (above), e.g., gala^xies with larger B/T do not have smaller Xre- Also, there is no 
correlation between Aije and Sersic index n, as shown in Figure 9 (below), e.g. galaxies with disk-like pseudo-bulges (with 
small n) do not have a tendency to be more rotationally supported. 

Good news in the models by Khochfar et al. [41] is also that it predicts a fairly large range of B/Ts for the fast 
rotators. However, the used B/T range of 0.5-0.9 is almost completely outside the range obtained by the multi-component 
decompositions for the SOs (see Fig. 8). In the literature possible mechanisms have been discussed which might lower the B/T 
in SAMS: (a) increasing the energy injected by supernovae ([2]), or (b) taking into account mass loss from already existing 
stars, which can lower B/T by a factor of 2-3 ([56]). However, these explanations have been criticized by Scannapieco et al. 
[65]. The first because during the mass loss the bulge density is also reduced, and the density of the disk increases only in 
the central parts of the disk. And the second, because the suggested efficiency of SN feedback is too high compared to the 
cosmological models, e.g., by Guo et al. [61]. In the current models by Khochfar et al. lower B/T values for fast rotators 
are obtained by significantly decreasing the frequency of minor mergers, but that might not be consistent with the observed 
merger frequency in the Universe. 



^ The 2-component decompositions arc thoroughly discussed in the recent review by Graham, to appear in the next edition of "Planets, 
Stars, and Stellar Systems", Vol 6, 2011. 
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5.4 Low B/T in conjunction with tidal models (e.g. Bekki & Couch [11]) 

In the stripping models the B/T problem has generally been the opposite: it has been difficult to make such massive bulges as 
observed in galaxies. In Bekki & Couch [11] the problem was solved by assuming that the tidal effects increase the efficiency 
of star formation and might also add some fresh material into the galaxies during the encounter process. Here we discuss how 
well the B/T ratios obtained in NIRSOS are compatible with these models. 

Compared to Khochfar et al. [41], in Bekki & Couch [11] a lower B/T range of 0.3-0.5 was used for the SOs. This is well 
within the range of the measured B/T values in Figure 8, but still lacking the lower B/T tail. In the simulations by Bekki & 
Couch, B/T is measured kinematically from the simulations by linking each simulated star either to a bulge or a disk. Based 
on the comparison by Scannapieco et al. [65], such kinematic decompositions give significantly higher B/T values than the 
decomposition of the flux distribution in galaxies, done in a similar manner as in observations. A similar conclusion was made 
also by Laurikainen & Salo [85]. Taking this into account the B/T range in the models by Bekki & Couch might be even more 
closer with the observed B/T range. 



6 PROPERTIES OF BULGES AND DISKS, BASED ON THE DECOMPOSITIONS 
6.1 Properties of bulges and disks: SOs vs. spirals 

The photometric properties of bulges and disks were compared between SOs and spirals by Laurikainen et al. [53], without 
including any kinematic information. The parameters were corrected for the internal dust in galaxies. The final sample, with 
reliable decompositions made, consists of 340 galaxies. Naively, strong correlations between the parameters of the bulge and 
the disk are expected if the disks were formed first, and the bulges emerged from the disc by secular evolution. On the other 
hand, in hierarchical clustering, followed by galaxy mergers, the properties of bulges were established already during the last 
major merger event ([39]), and presumably changed very little after that. Also, if the formative processes of bulges in SOs 
were similar to those of the elliptical galaxies, a continuity in the properties of bulges and ellipticals is expected. 

It was found that the bulges of SOs are different from the elliptical galaxies, and similar to the bulges in spirals, manifested 
in the r^ff (bulge) vs. Mb^ige diagram (Fig. 10): bulges both for the SOs and spirals appear below the line for the elliptical 
galaxies, indicated by the dashed lin^. It was shown by Graham & Worley ([35], their Fig. 11) that this line nicely follows 
the observed distribution of the ellipticals (see [34]), covering a large magnitude range, and appearing in various types of 
environments. For spirals the discontinuity of bulges and ellipticals was shown also by Gadotti [33]. The fact that the bulges 
of SOs are not similar to the elliptical galaxies is manifested also in the Kormendy relation and in the photometric plane where 
only the brightest bulges of SOs follow the location of the Coma cluster ellipticals (see Figs. 2 and 3 in [53]). 

The rf.f / (bulge) vs. Mbuige in Figure 10 extends to the regime of very small bulges, corresponding to the effective radii of 
the dwarf early-type galaxies, dEs. In Figure 10 the dEs would appear in the flat part of the dashed line (see [34]). However, 
it is worth noticing that the small bulges in SOs have several magnitudes higher surface brightnesses than dEs, and therefore 
do not fit to the scenario in which dEs were formed from Sc type spirals by loosing their disks in galaxy harassment in dense 
cluster potentials. This is consistent with the recent study by Kormendy & Bender [88], who suggested that the red and dead 
dwarf galaxies are rather transferred from later type Scd-Im galaxies. 

Not only bulges, but also the disks in SOs are similar to those in spirals. It was shown by Laurikainen et al. [53], that 
both for the SOs and spirals the scale length of the disk (hn) varies between 1-10 kpc, and the central surface brightness {fio) 
between 16-20 magnitudes. Also, fio has no correlation with the mass of the disk {Mdisk)- It was concluded that does not 
follow the Freeman law ([32]) pertaining a constant central surface brightness (see Fig. 9 in [53]). The similarity in the disc 
sizes between SOs and spirals is consistent with the tidal models by Bekki & Couch [11] for the origin of SOs. They stated that 
"it is unlikely that SOs have systematically larger disk masses in comparison with spirals." 

For the SOs the photometric properties of bulges and disks in NIRSOS are found to be coupled, in a similar manner as for 
spirals [53]. This is manifested in correlations between their scale parameters, hn vs. r-g// (bulge), and their masses, Mbuige vs. 
Mdisk, as shown in Figure 11. Also, for the bulges having Mbuige ~ -20 mag in the Ks-hand the difference in Mdisk between 
the early (SO, SO/a, Sa) and late type (Sab and later) galaxies is nearly one magnitude. Based on the stellar population 
models by Bedregal, Aragon-Salamanca & Merrifield [10] this difference could be explained by quiescent star formation if star 
formation of the disc in spiral galaxies was truncated 1-6 Gyr ago. 

These results are consistent with the picture in which the discs were formed first and the bulges emerged from the disc 
material afterwords, without invoking any merger events, which might destroy the observed correlations between the bulges 
and disks. However, our results do not unambiguously rule out the other scenarios. For example, in minor mergers the satellite 
galaxies increase the bulge mass in one hand, and at the same time the disks stars are spilled to a larger radius, which might 
also lead to a correlation between the scale parameters of the bulge and the disk. 

^ the dashed line is not for the Coma cluster ellipticals, as erroneously claimed in Laurikainen et al. [53] 
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6.2 Properties of bulges in the Hubble sequence 

Figure 12 shows the mean total galaxy brightnesses in the i^s-band and the parameters of the bulge in the Hubble sequence 
[53]. The bulge parameters are compared with those by Graham & Worley [35]. Notice that the comparison is meaningful 
only between T=0-5 where enough common galaxies appear in both studies. 

The following things can be noticed: (1) in the near-IR, tracing the old stellar population of galaxies, the SOs in NIRSOS 
are not less luminous than the spirals in OSUBSGS (both samples are magnitude- limited, selected in a similar manner). A 
similar result in the near-IR was obtained also by Burnstein ot al. ([89]). This is a puzzle if SOs were formed from spirals Sc or 
earlier, simply by loosing their gas in the stripping mechanisms, followed by quenching of star formation. The luminosities of 
SOs could be explained if they have gained more mass after the SO characteristics were settled down, for example by accreting 
a significant number of satellite galaxies during their lifetime. However, that is not supported by our analysis of the Qg 
distribution of bars in SOs, at least if the accreting satellites were gas rich. A more likely solution to this puzzle is that at 
higher redshifts some fraction of the progenitors of SOs were more luminous than the local spirals. This is evidenced by Geach 
et al. [90] showing that LIRGS at z=0.5 could account for the star formation needed to explain the bright end of the SO galaxy 
luminosity function (see the discussion in [53]). 

(2) Veff/hR decreases slightly from SOs towards spirals. This parameter is smaller than obtained in the previous studies, 
mainly because we use multi-component decompositions which presumes the erroneous additions of bars to the bulges. For all 
Hubble types the observed mean reff/hn is smaller than predicted by the current merger models (e.g. [58]). (3) The Sersic 
index n is very similar for SO, SO/a and Sa galaxies, and slightly decreases towards the later types. Such small values of Sersic 
index (n ^ 2) are not predicted in the current merger models (e.g. [66], [59]). (4) B/T was already discussed in Section 5. 
Figure 12 additionally shows that the mean B/T is very similar for SO, SO/a and Sa galaxies, and only after that decreases 
towards the later types. It was also shown by Laurikainen et al. [51] that B/T is lower for barred than non-barred galaxies, 
and also lower for galaxies having flat-top DG bars compared to the more simple SG bars. 

Overall, the above properties of bulges are not well compatible with the current merger models. 



7 THE OUTLINED PICTURE 

7.1 SOs transformed from spirals: the Hubble sequence revisited 

The Hubble sequence is revisited, following the early suggestion by Spitzer & Baade [71] and van den Bergh [72]: namely the 
SOs are suggested to be spread throughout the Hubble sequence in parallel tuning forks as spirals (SOa, SOb, SOc etc.) ([53], 
[54]). Originally this scenario was arisen by van den Bergh because he found that the entire SO class was offset towards 
fainter optical magnitudes from the luminosity distribution of Sa galaxies. He then assumed that there must exist anemic SOs 
whose surface brightness distributions are similar to the faint surface brightness levels of Sb and Sc galaxies. However, at that 
time no such galaxies were found. Our NIRSOS/OSUBSGS analysis has shown that indeed, the disks in SOs have a similar 
range of central surface brightnesses and scale lengths of the disk as the disks of Sa-Sc type spirals. And in particular, our 
multi-component decomposition approach has shown that the bulges in many SOs have B/T ratios as small as typically found 
in Sc typo spirals [53]. A similar view to the Hubble sequence was taken recently also by Cappellari ot al. [22], based on new 
kinematic observations, showing that bulges in SOs are fast rotating oblate systems. During the referee process of this work a 
paper by Kormendy & Bender [88] appeared in astro-ph, in which a similar view was also suggested. 

Our results are consistent with the idea according to which SOs are mainly former spirals transformed into SOs by loosing 
their gas in some stripping mechanisms. The following results of our NIRSOS/OSUBSGS analysis are consistent with this 
picture: (a) SOs cover the full observed S/T-range of spirals, (b) The photometric properties of bulges and disks appeared to 
be similar in SOs and spirals. And most importantly, the bulges in SOs have no continuity with the elliptical galaxies in the 
Te// (bulge) vs. Mbuige diagram, (c) The relative masses and scale parameters of the bulges and disks in SOs are coupled, in 
a similar manner as in spirals, which hint to a common origin of the bulges and disks in these galaxies. What needs to be 
explained in this scenario is why the SOs in NIRSOS are not less luminous than their assumed spiral progenitors. 

7.2 Solution for the bulge problem in the stripping scenario? 

In the stripping scenario, a problem of larger bulges in SOs, compared to bulges in spirals, was discussed. It was concluded, 
based on the stellar population models by Bedregal, Aragon-Salamanca & Merrifield [10], that spirals with bulges brighter 
than -20 magnitudes in the Jf-band are massive enough to be converted into bulges of SOs, simply by slow internal galax;y 
evolution. However, that cannot explain why SOs are not dimmer than the spirals, as expected if they lost some of their gas 
in the stripping process. This implies that some other mechanism is also needed to complete the mass accretion in the bulges 
of SOs. Possible alternatives are that the progenitors of SOs were more gas rich than the present day spirals, which would 
explain also the similar total galaxy brightnesses in most SOs, compared to the galaxy brightnesses in spirals, in the near-IR. 
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Or some mechanism other than gas accretion is taking place. One suggested mechanism is the so called potential-density 
phase shift mechanism suggested by Zhang [74], allowing stellar mass transfer even after the gas depletion has stopped. Or, 
as suggested by Bekki & Couch [11], tidal effects in galaxy groups can increase star formation efficiency in the central regions. 
Minor mergers can also accret material to the main galaxies, but that process needs to be mild enough for not to increase the 
Sersic index over the fairly low observed values, at least in the current models. Also, if our conclusion of the Qg distribution 
for the SOs is correct, the accreting satellites need to be fairly gas poor. 

7.3 Bars and lenses: evidence of secular evolution? 

A test was made for the simulation models by Bournaud & Combes [14] , in which bars are expected to have multiple episodes 
during the Hubble time. Based on the observed bar strength (Qg) distribution we concluded that the SOs have not accreted 
gas for a long time. Spirals show a larger number of strong bars, but the observed Qg distribution even for them resembles 
the non-accretion models. This implies that most probably bars are fairly robust in all Hubble types, without experiencing 
multiple bar episodes in the Hubble time. However, inspite of that bars may evolve significantly over time, as manifested in 
the morphology of the old stellar population of bars in the near-IR. Namely bars in SOs have ansae morphology and double- 
peaked Fourier density profiles much more frequently than bars in spiral galaxies. According to the dynamical models such 
features can be associated to evolved strong bars. Bars may also dissolve into lenses, as first suggested by Kormendy [43], 
which is consistent with the observed lower bar fraction, and a larger fraction of lenses among the SOs. However, based on the 
morphology of lenses, not all lenses can be dissolved former bars. 

A new lens type, barlens, was introduced, appearing in strong bars, most probably forming part of the bar itself. We 
assume that it might be related to the re-distribution of matter in galaxies. Also, X-shaped bars in nearly face-on galaxies 
were detected, which cannot be simply vertically thickened bars by buckling effects, and therefore needs to be explained in the 
framework of the theoretical models. Multiple lenses were detected even in 25% of the SOs, including barred and non-barred 
galaxies, which are not yet well understood in the current paradigm of galaxy formation and evolution. We anticipate that 
secular evolution might play an important role in the formative processes of lenses. 
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Figure 1. An example of the Atlas image [54]. The upper panel shows the whole galaxy using the full magnitude range, for which the 
scale is indicated in the right. The scales in x- and y-axis are in arcseconds. The small panels show the same image, but selecting the 
radial and magnitude scales in a different manner: the two upper panels use the full radial scale, whereas in the lower panel only the 
central region of the galaxy is shown. The morphological classification is from the Atlas paper. The arrow indicates the barlens. The 
graphics show the radial profiles of the surfax;e brightness (upper panel), ellipticity (second panel), position angle (third upper panel), 
and parameter b4 (lowermost panel) indicating deviations from perfect ellipticity. The logarithmic radial scale is in arcseconds. In these 
plots the vertical dashed lines indicate the bar radius, and the full line is the radius where the surface brightness is 20 mag arcseconds"^. 





Figure 2. This figure is compiled of two figures that originally appeared in Laurikainen et al. [54]. The images are flux calibrated and are 
shown in a logarithmic scale, covering the full magnitude range of the images. In the upper row the original images are shown, whereas 
in the lower row the bar and bulge models from the decomposition are subtracted from the original image of NGC 2983. These images 
highlight the main type of bar morphologies in SO galaxies. 
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NGC 1 543 (l,nl) central parts (enlarged by 1 0) 




Figure 3. Two panels from Figure 15 in Laurikainen et al. [54] are shown. In the left panel the whole galaxy is shown using the full 
magnitude range of the image. In this image the outer lens and the inner lens surrounding the main bar are clearly visible. In the right 
panel only the central regions of the galaxy is shown (enlarged by a factor of 10) , demonstrating the morphology of the nuclear bar, and 
the nuclear lens extending to the radius of the nuclear bar. 



NIRSOS review 13 



NGC 3941 NGC 4245 




r(arcsec) r(arcsec) 



Figure 4. Radial profiles of the Fourier amplitudes of density: left panel shows a simple bar, whereas right panel shows a two-component 

bar: both lower and higher Fourier modes are significant in these two components. The thin lines show the m=2, 4, 6, 8 Fourier amplitudes, 
normalized to m=0. The thick grey line shows the force ratio QT(r) = \FT{r, (j>)\max/ < \FR{r, (f))\ > (Qg indicates the maximum of QT 
), and the vertical line indicates the length of the bar. The crosses show individual Gaussian fits to the m=2 Fourier amplitude profile. 
The single peaked bar represents a thin classical bar. The double peaked amplitude profile is associated to a bar having a shorter and 
vertically thicker inner part, and a longer and vertically thinner outer part. This figure is a modified version of Figure 8 in Laurikainen 
et al. [51]. 
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Figure 5. Left: the observed Qg distribution for spirals in the OSUBSGS (cover the Hubble types up to T=9), tciken from Block et al. 
[80], compared with the simulation models by Bournaud & Combes [14]. The solid line shows the observed values from Block et al., the 
dotted line is the model with no gas a<;cretion, and the dashed line shows the model that doubles the total galaxy mass in 10^" years. 
Right: the observed Qg distributions for the SOs and early- type spirals in NIRSOS, talcen from Buta et al. [86]. Notice that Qf, in Block 
et al. [80] in the left panel has the same meaning as Qg in Buta et al. [86] in the right panel. 
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Figure 6. Examples of non-barred lens galaxies. Left panels show the images, whereas right panels show the 2-dimensional surface 
brightness profiles, in which over-plotted are the fitted functions to the structural decompositions, as well as the total model. White dots 
are for the image pixel values, and the black lines/dots are for the total model. Models for the bulges are shown by dark grey lines, and 
the models for the lenses by light grey and black dots in the inner regions. For example, in NGC 524 the two curved inner components 
are the nuclear lens (nl) and lens (1), and the outer exponential function fits the outer lens (L). This figure was originally published by 
Laurikainen et al. [54]. 
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Figure 7. Upper' row: surface brightness ijrofilcs of synthetic images in which noise is added to mimic a real image. In a) the image 
has only a bulge and a disk, whereas in b) and c) a small bar is added on top of the disk. Black dots show the pixel values of the 
synthetic images, and the other symbols show the functions fitted to the 2-dimensional images, as explained in the text. Lower row: 
the residual images for the above decompositions are shown. They are obtained by subtracting the total decomposition model from the 
original synthetic image. The figure is a modified version of Figure 2 in Laurikainen et al. [49]. 
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Figure 8. Histogram of the B/T flux ratios in Ks-haiid. corrected for Galactic and internal dust extinction. The corrections are different 
for the bulge and the disk. The histograms for the SOs (full line), and the SO/a galaxies (dotted line) are based on the NIRSOS sample. 
For spirals the lines for Sab-Sbc (full white line) and later than Sc (dashed white line) are drawn using the OSUBSGS sample. This 
figure was originally published by Laurikainen et al. [53]. 
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Figure 9. The kinematic parameter A_r vs. B/T {upper panel), and vs. Sersic n {lower panel). Parameter Xr^, taJcen from Emsellem et 
al. [29] can be used as a discriminator between oblate and triaxial nature of the galaxies, within the apertures used in the measurements. 
In this plot Ah obtained within one effective radius of the galaxy is used. The n-parameters and the dust corrected B/T values are taken 
from Laurikainen et al. [53] . 
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Figure 10. The eflfective radius of the bulge (Vf-ff) plotted against the absolute brightness of the bulge (M^(imige)). Shown separately 
are: SOs (filled symbols), Sa-Sbc spirals (crosses), and Sc and later types (squares). The dashed line shows the location of the elliptical 
galaxies, taken from Graham & Worley (2008). The continuous lines show the mean values of the data points in one magnitude bins: the 
thick line is for the SOs and the thin line for the Sab-Sbc spirals. This figure is part of Figure 8a in Laurikainen et al. [53]. 
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Figure 11. Left: the absolute brightness of the disk {M'^(disk)) is plotted against the absolute brightness of the bulge (M^Qmlge)). 
(a) SOs (red filled circles), SO/a+Sa (open circles), Sab-Sbc (crosses) and Sc and later types (squares). Solid, dash-dotted and dashed 
lines show linear fits for the SOs, SO/a+Sa and Sab-Sbc galaxies, respectively. All the correlations are statistically significant. Right: hn 
(disk) vs. r^ff (bulge) for the galaxies in NIRSOS, divided to two galaxy brightness bins in i^a-band. This figure is compiled from two 
figures originally published by Laurikainen et al. [52] , [53] . 
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Figure 12. As a function of Hubble type are shown: (a) the dust corrected B/T, (b) Scrsic index n, (c) the ratio of the effective radius 
of the bulge (r^ff) divided by the scale length of the disk {h°), and (d) the total absolute galaxy magnitude corrected for Galactic 
extinction (M^^^.). The small symbols show the measurements for the individual galaxies (with a small random horizontal offset), filled 
circles show the median values in each Hubble type bin, and large open circles are the median values from Graham & Worley ([35]; 
their value for SOs refers to types -3, -2, -1, and is here drawn at T=-2). The error bars are taJjen to be the sample variances in each 
bin divided by square root of the number of galaxies in each bin. Dark symbols are for the galaxies in NIRSOS and light symbols for the 
galaxies in OSUBSGS. In a slightly different form the same figure was originally published by Laurikainen et al. [53]. 
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